\(\int \frac {(a d e+(c d^2+a e^2) x+c d e x^2)^{3/2}}{(d+e x)^4} \, dx\) [1926]

Optimal result

Integrand size = 37, antiderivative size = 169 \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{(d+e x)^4} \, dx=-\frac {2 c d \sqrt {a d e+\left (c d^2+a e^2\right ) x+c d e x^2}}{e^2 (d+e x)}-\frac {2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{3 e (d+e x)^3}+\frac {c^{3/2} d^{3/2} \text {arctanh}\left (\frac {c d^2+a e^2+2 c d e x}{2 \sqrt {c} \sqrt {d} \sqrt {e} \sqrt {a d e+\left (c d^2+a e^2\right ) x+c d e x^2}}\right )}{e^{5/2}} \]

[Out]

-2/3*(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^(3/2)/e/(e*x+d)^3+c^(3/2)*d^(3/2)*arctanh(1/2*(2*c*d*e*x+a*e^2+c*d^2)/c
^(1/2)/d^(1/2)/e^(1/2)/(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^(1/2))/e^(5/2)-2*c*d*(a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2
)^(1/2)/e^2/(e*x+d)

Rubi [A] (verified)

Time = 0.06 (sec) , antiderivative size = 169, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.081, Rules used = {676, 635, 212} \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{(d+e x)^4} \, dx=\frac {c^{3/2} d^{3/2} \text {arctanh}\left (\frac {a e^2+c d^2+2 c d e x}{2 \sqrt {c} \sqrt {d} \sqrt {e} \sqrt {x \left (a e^2+c d^2\right )+a d e+c d e x^2}}\right )}{e^{5/2}}-\frac {2 c d \sqrt {x \left (a e^2+c d^2\right )+a d e+c d e x^2}}{e^2 (d+e x)}-\frac {2 \left (x \left (a e^2+c d^2\right )+a d e+c d e x^2\right )^{3/2}}{3 e (d+e x)^3} \]

[In]

Int[(a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2)^(3/2)/(d + e*x)^4,x]

[Out]

(-2*c*d*Sqrt[a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2])/(e^2*(d + e*x)) - (2*(a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x
^2)^(3/2))/(3*e*(d + e*x)^3) + (c^(3/2)*d^(3/2)*ArcTanh[(c*d^2 + a*e^2 + 2*c*d*e*x)/(2*Sqrt[c]*Sqrt[d]*Sqrt[e]
*Sqrt[a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2])])/e^(5/2)

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 635

Int[1/Sqrt[(a_) + (b_.)*(x_) + (c_.)*(x_)^2], x_Symbol] :> Dist[2, Subst[Int[1/(4*c - x^2), x], x, (b + 2*c*x)
/Sqrt[a + b*x + c*x^2]], x] /; FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 676

Int[((d_.) + (e_.)*(x_))^(m_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(d + e*x)^(m + 1)*((
a + b*x + c*x^2)^p/(e*(m + p + 1))), x] - Dist[c*(p/(e^2*(m + p + 1))), Int[(d + e*x)^(m + 2)*(a + b*x + c*x^2
)^(p - 1), x], x] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && EqQ[c*d^2 - b*d*e + a*e^2, 0] && GtQ[
p, 0] && (LtQ[m, -2] || EqQ[m + 2*p + 1, 0]) && NeQ[m + p + 1, 0] && IntegerQ[2*p]

Rubi steps \begin{align*} \text {integral}& = -\frac {2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{3 e (d+e x)^3}+\frac {(c d) \int \frac {\sqrt {a d e+\left (c d^2+a e^2\right ) x+c d e x^2}}{(d+e x)^2} \, dx}{e} \\ & = -\frac {2 c d \sqrt {a d e+\left (c d^2+a e^2\right ) x+c d e x^2}}{e^2 (d+e x)}-\frac {2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{3 e (d+e x)^3}+\frac {\left (c^2 d^2\right ) \int \frac {1}{\sqrt {a d e+\left (c d^2+a e^2\right ) x+c d e x^2}} \, dx}{e^2} \\ & = -\frac {2 c d \sqrt {a d e+\left (c d^2+a e^2\right ) x+c d e x^2}}{e^2 (d+e x)}-\frac {2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{3 e (d+e x)^3}+\frac {\left (2 c^2 d^2\right ) \text {Subst}\left (\int \frac {1}{4 c d e-x^2} \, dx,x,\frac {c d^2+a e^2+2 c d e x}{\sqrt {a d e+\left (c d^2+a e^2\right ) x+c d e x^2}}\right )}{e^2} \\ & = -\frac {2 c d \sqrt {a d e+\left (c d^2+a e^2\right ) x+c d e x^2}}{e^2 (d+e x)}-\frac {2 \left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{3 e (d+e x)^3}+\frac {c^{3/2} d^{3/2} \tanh ^{-1}\left (\frac {c d^2+a e^2+2 c d e x}{2 \sqrt {c} \sqrt {d} \sqrt {e} \sqrt {a d e+\left (c d^2+a e^2\right ) x+c d e x^2}}\right )}{e^{5/2}} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.29 (sec) , antiderivative size = 130, normalized size of antiderivative = 0.77 \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{(d+e x)^4} \, dx=\frac {2 \sqrt {(a e+c d x) (d+e x)} \left (-\frac {\sqrt {e} \left (a e^2+c d (3 d+4 e x)\right )}{(d+e x)^2}+\frac {3 c^{3/2} d^{3/2} \text {arctanh}\left (\frac {\sqrt {c} \sqrt {d} \sqrt {d+e x}}{\sqrt {e} \sqrt {a e+c d x}}\right )}{\sqrt {a e+c d x} \sqrt {d+e x}}\right )}{3 e^{5/2}} \]

[In]

Integrate[(a*d*e + (c*d^2 + a*e^2)*x + c*d*e*x^2)^(3/2)/(d + e*x)^4,x]

[Out]

(2*Sqrt[(a*e + c*d*x)*(d + e*x)]*(-((Sqrt[e]*(a*e^2 + c*d*(3*d + 4*e*x)))/(d + e*x)^2) + (3*c^(3/2)*d^(3/2)*Ar
cTanh[(Sqrt[c]*Sqrt[d]*Sqrt[d + e*x])/(Sqrt[e]*Sqrt[a*e + c*d*x])])/(Sqrt[a*e + c*d*x]*Sqrt[d + e*x])))/(3*e^(
5/2))

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(472\) vs. \(2(147)=294\).

Time = 2.94 (sec) , antiderivative size = 473, normalized size of antiderivative = 2.80

method	result	size
default	\(\frac {-\frac {2 \left (c d e \left (x +\frac {d}{e}\right )^{2}+\left (e^{2} a -c \,d^{2}\right ) \left (x +\frac {d}{e}\right )\right )^{\frac {5}{2}}}{3 \left (e^{2} a -c \,d^{2}\right ) \left (x +\frac {d}{e}\right )^{4}}+\frac {2 c d e \left (-\frac {2 \left (c d e \left (x +\frac {d}{e}\right )^{2}+\left (e^{2} a -c \,d^{2}\right ) \left (x +\frac {d}{e}\right )\right )^{\frac {5}{2}}}{\left (e^{2} a -c \,d^{2}\right ) \left (x +\frac {d}{e}\right )^{3}}+\frac {4 c d e \left (\frac {2 \left (c d e \left (x +\frac {d}{e}\right )^{2}+\left (e^{2} a -c \,d^{2}\right ) \left (x +\frac {d}{e}\right )\right )^{\frac {5}{2}}}{\left (e^{2} a -c \,d^{2}\right ) \left (x +\frac {d}{e}\right )^{2}}-\frac {6 c d e \left (\frac {\left (c d e \left (x +\frac {d}{e}\right )^{2}+\left (e^{2} a -c \,d^{2}\right ) \left (x +\frac {d}{e}\right )\right )^{\frac {3}{2}}}{3}+\frac {\left (e^{2} a -c \,d^{2}\right ) \left (\frac {\left (2 c d e \left (x +\frac {d}{e}\right )+e^{2} a -c \,d^{2}\right ) \sqrt {c d e \left (x +\frac {d}{e}\right )^{2}+\left (e^{2} a -c \,d^{2}\right ) \left (x +\frac {d}{e}\right )}}{4 c d e}-\frac {\left (e^{2} a -c \,d^{2}\right )^{2} \ln \left (\frac {\frac {e^{2} a}{2}-\frac {c \,d^{2}}{2}+c d e \left (x +\frac {d}{e}\right )}{\sqrt {c d e}}+\sqrt {c d e \left (x +\frac {d}{e}\right )^{2}+\left (e^{2} a -c \,d^{2}\right ) \left (x +\frac {d}{e}\right )}\right )}{8 c d e \sqrt {c d e}}\right )}{2}\right )}{e^{2} a -c \,d^{2}}\right )}{e^{2} a -c \,d^{2}}\right )}{3 \left (e^{2} a -c \,d^{2}\right )}}{e^{4}}\)	\(473\)

[In]

int((a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^(3/2)/(e*x+d)^4,x,method=_RETURNVERBOSE)

[Out]

1/e^4*(-2/3/(a*e^2-c*d^2)/(x+d/e)^4*(c*d*e*(x+d/e)^2+(a*e^2-c*d^2)*(x+d/e))^(5/2)+2/3*c*d*e/(a*e^2-c*d^2)*(-2/
(a*e^2-c*d^2)/(x+d/e)^3*(c*d*e*(x+d/e)^2+(a*e^2-c*d^2)*(x+d/e))^(5/2)+4*c*d*e/(a*e^2-c*d^2)*(2/(a*e^2-c*d^2)/(
x+d/e)^2*(c*d*e*(x+d/e)^2+(a*e^2-c*d^2)*(x+d/e))^(5/2)-6*c*d*e/(a*e^2-c*d^2)*(1/3*(c*d*e*(x+d/e)^2+(a*e^2-c*d^
2)*(x+d/e))^(3/2)+1/2*(a*e^2-c*d^2)*(1/4*(2*c*d*e*(x+d/e)+e^2*a-c*d^2)/c/d/e*(c*d*e*(x+d/e)^2+(a*e^2-c*d^2)*(x
+d/e))^(1/2)-1/8*(a*e^2-c*d^2)^2/c/d/e*ln((1/2*e^2*a-1/2*c*d^2+c*d*e*(x+d/e))/(c*d*e)^(1/2)+(c*d*e*(x+d/e)^2+(
a*e^2-c*d^2)*(x+d/e))^(1/2))/(c*d*e)^(1/2))))))

Fricas [A] (verification not implemented)

none

Time = 0.47 (sec) , antiderivative size = 426, normalized size of antiderivative = 2.52 \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{(d+e x)^4} \, dx=\left [\frac {3 \, {\left (c d e^{2} x^{2} + 2 \, c d^{2} e x + c d^{3}\right )} \sqrt {\frac {c d}{e}} \log \left (8 \, c^{2} d^{2} e^{2} x^{2} + c^{2} d^{4} + 6 \, a c d^{2} e^{2} + a^{2} e^{4} + 4 \, {\left (2 \, c d e^{2} x + c d^{2} e + a e^{3}\right )} \sqrt {c d e x^{2} + a d e + {\left (c d^{2} + a e^{2}\right )} x} \sqrt {\frac {c d}{e}} + 8 \, {\left (c^{2} d^{3} e + a c d e^{3}\right )} x\right ) - 4 \, \sqrt {c d e x^{2} + a d e + {\left (c d^{2} + a e^{2}\right )} x} {\left (4 \, c d e x + 3 \, c d^{2} + a e^{2}\right )}}{6 \, {\left (e^{4} x^{2} + 2 \, d e^{3} x + d^{2} e^{2}\right )}}, -\frac {3 \, {\left (c d e^{2} x^{2} + 2 \, c d^{2} e x + c d^{3}\right )} \sqrt {-\frac {c d}{e}} \arctan \left (\frac {\sqrt {c d e x^{2} + a d e + {\left (c d^{2} + a e^{2}\right )} x} {\left (2 \, c d e x + c d^{2} + a e^{2}\right )} \sqrt {-\frac {c d}{e}}}{2 \, {\left (c^{2} d^{2} e x^{2} + a c d^{2} e + {\left (c^{2} d^{3} + a c d e^{2}\right )} x\right )}}\right ) + 2 \, \sqrt {c d e x^{2} + a d e + {\left (c d^{2} + a e^{2}\right )} x} {\left (4 \, c d e x + 3 \, c d^{2} + a e^{2}\right )}}{3 \, {\left (e^{4} x^{2} + 2 \, d e^{3} x + d^{2} e^{2}\right )}}\right ] \]

[In]

integrate((a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^(3/2)/(e*x+d)^4,x, algorithm="fricas")

[Out]

[1/6*(3*(c*d*e^2*x^2 + 2*c*d^2*e*x + c*d^3)*sqrt(c*d/e)*log(8*c^2*d^2*e^2*x^2 + c^2*d^4 + 6*a*c*d^2*e^2 + a^2*
e^4 + 4*(2*c*d*e^2*x + c*d^2*e + a*e^3)*sqrt(c*d*e*x^2 + a*d*e + (c*d^2 + a*e^2)*x)*sqrt(c*d/e) + 8*(c^2*d^3*e
+ a*c*d*e^3)*x) - 4*sqrt(c*d*e*x^2 + a*d*e + (c*d^2 + a*e^2)*x)*(4*c*d*e*x + 3*c*d^2 + a*e^2))/(e^4*x^2 + 2*d
*e^3*x + d^2*e^2), -1/3*(3*(c*d*e^2*x^2 + 2*c*d^2*e*x + c*d^3)*sqrt(-c*d/e)*arctan(1/2*sqrt(c*d*e*x^2 + a*d*e
+ (c*d^2 + a*e^2)*x)*(2*c*d*e*x + c*d^2 + a*e^2)*sqrt(-c*d/e)/(c^2*d^2*e*x^2 + a*c*d^2*e + (c^2*d^3 + a*c*d*e^
2)*x)) + 2*sqrt(c*d*e*x^2 + a*d*e + (c*d^2 + a*e^2)*x)*(4*c*d*e*x + 3*c*d^2 + a*e^2))/(e^4*x^2 + 2*d*e^3*x + d
^2*e^2)]

Sympy [F(-1)]

Timed out. \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{(d+e x)^4} \, dx=\text {Timed out} \]

[In]

integrate((a*d*e+(a*e**2+c*d**2)*x+c*d*e*x**2)**(3/2)/(e*x+d)**4,x)

[Out]

Timed out

Maxima [F(-2)]

Exception generated. \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{(d+e x)^4} \, dx=\text {Exception raised: ValueError} \]

[In]

integrate((a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^(3/2)/(e*x+d)^4,x, algorithm="maxima")

[Out]

Exception raised: ValueError >> Computation failed since Maxima requested additional constraints; using the 'a
ssume' command before evaluation *may* help (example of legal syntax is 'assume(e*(a*e^2-c*d^2)>0)', see `assu
me?` for mor

Giac [F(-2)]

Exception generated. \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{(d+e x)^4} \, dx=\text {Exception raised: TypeError} \]

[In]

integrate((a*d*e+(a*e^2+c*d^2)*x+c*d*e*x^2)^(3/2)/(e*x+d)^4,x, algorithm="giac")

[Out]

Exception raised: TypeError >> an error occurred running a Giac command:INPUT:sage2:=int(sage0,sageVARx):;OUTP
UT:Unable to divide, perhaps due to rounding error%%%{%%{[%%%{1,[0,0,10]%%%},0]:[1,0,%%%{-1,[1,1,1]%%%}]%%},[4
,4]%%%}+%%%

Mupad [F(-1)]

Timed out. \[ \int \frac {\left (a d e+\left (c d^2+a e^2\right ) x+c d e x^2\right )^{3/2}}{(d+e x)^4} \, dx=\int \frac {{\left (c\,d\,e\,x^2+\left (c\,d^2+a\,e^2\right )\,x+a\,d\,e\right )}^{3/2}}{{\left (d+e\,x\right )}^4} \,d x \]

[In]

int((x*(a*e^2 + c*d^2) + a*d*e + c*d*e*x^2)^(3/2)/(d + e*x)^4,x)

[Out]

int((x*(a*e^2 + c*d^2) + a*d*e + c*d*e*x^2)^(3/2)/(d + e*x)^4, x)